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Abstract. Using the United Nations COMTRADE database [T] we construct the Google matrix G of 
multiproduct world trade between the UN countries and analyze the properties of trade flows on this 
network for years 1962 - 2010. This construction, based on Markov chains, treats all countries on equal 
democratic grounds independently of their richness and at the same time it considers the contributions of 
trade products proportionally to their trade volume. We consider the trade with 61 products for up to 227 
countries. The obtained results show that the trade contribution of products is asymmetric: some of them 
are export oriented while others are import oriented even if the ranking by their trade volume is symmetric 
in respect to export and import after averaging over all world countries. The construction of the Google 
matrix allows to investigate the sensitivity of trade balance in respect to price variations of products, e.g. 
petroleum and gas, taking into account the world connectivity of trade links. The trade balance based on 
PageRank and GheiRank probabilities highlights the leading role of Ghina and other BRIGS countries in 
the world trade in recent years. We also show that the eigenstates of G with large eigenvalues select specific 
trade communities. 

PACS. 89.75.Fb Structures and organization in complex systems - 89.65.Gh Econophysics - 89.75.He 
Networks and genealogical trees - 89.20.Hh World Wide Web, Internet 


1 Introduction 

According to the data of UN COMTRADE [Tj and the in¬ 
ternational trade statistics 2014 of the World Trade Orga¬ 
nization (WTO) [2] the international world trade between 
world countries demonstrates a spectacular growth with 
an increasing trade volume and number of trade products. 
It is well clear that the world trade plays the fundamental 
role in the development of world economy [3]. According 
to the WTO Chief Statistician Hubert Escaith “In recent 
years we have seen growing demand for data on the world 
economy and on international trade in particular. This 
demand has grown in particular since the 2008-09 crisis, 
whose depth and breadth surprised many experts” la.In 
global the data of the world trade exchange can be viewed 
as a large multi-functional directed World Trade Network 
(WTN) which provides important information about mul¬ 
tiproduct commercial flows between countries for a given 
year. At present the COMTRADE database contains data 
for Nc = 227 UN countries with up to Np ^ 10^ trade 
products. Thus the whole matrix of these directed trade 
flows has a rather large size N = NpNc r\j 10^. A usual 
approach is to consider the export and import volumes, 
expressed in US dollars (USD). An example of the world 
map of countries characterized by their import and export 
trade volume for year 2008 is shown in Eig. However, 
such an approach gives only an approximate description of 


trade where hidden links and interactions between certain 
countries and products are not taken into account since 
only a country global import or export are considered. 
Thus the statistical analysis of these multiproduct trade 
data requires a utilization of more advanced mathematical 
and numerical methods. 

In fact, in the last decade, modern societies developed 
enormous communication and social networks including 
the World Wide Web (WWW), Wikipedia, Twitter etc. 
(see e.g. 0 ). A necessity of information retrieval from 
such networks led to a development of efficient algorithms 
for information analysis on such networks appeared in 
computer science. One of the most spectacular tools is 
the PageRank algorithm developed by Brin and Page in 
1998 [6], which became a mathematical foundation of the 
Google search engine (see e.g. 0 )- This algorithm is based 
on the concept of Markov chains and a construction of the 
Google matrix G of Markov transitions between network 
nodes. The right eigenvector of this matrix G, known as 
PageRank vector, allows to rank all nodes according to 
their importance and influence on the network. The stud¬ 
ies of various directed networks showed that it is useful 
to analyze also the matrix G* constructed for the same 
network but with an inverted direction of links EMg. The 
PageRank vector of G* is known as the GheiRank vec¬ 
tor. The spectral properties of Goggle matrix for various 
networks are described in m- 
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Fig. 1. World map of countries with color showing country 
import (top panel) and export (bottom panel) trade volume 
expressed in millions of USD given by numbers of the color 
bars. The data are shown for year 2008 with Nc = 227 countries 
for trade in all Np = 61 products (from UN COMTRADE [T]). 
Names of countries can be hnd at [4]. 


The approach of Google matrix to the analysis of WTN 
was started in m- The striking feature of this approach 
is that it treats all UN countries on equal democratic 
grounds, independently of richness of a given country, in 
agreement with the principles of UN where all countries 
are equal. This property of G matrix is based on the prop¬ 
erty of Markov chains where the total probability is con¬ 
served to be unity since the sum of elements for each col¬ 
umn of G is equal to unity. Even if in this approach all 
countries are treated on equal grounds still the PageR- 
ank and CheiRank analysis recover about 75% of indus¬ 
trially developed countries of G20. However, now these 
countries appear at the top ranking positions not due to 
their richness but due to the efficiency of their trade net¬ 
work. Another important aspect found in m is that both 
PageRank and CheiRank vectors appear very naturally in 
the WTN corresponding to import and export flows. 

In this work we extend the Google matrix analysis for 
the multiproduct WTN obtained from COMTRADE [T] 
with up to A/p = 61 trade products for up to Nc = 227 
countries. The global G matrix of such trade flows has a 
size up to = NpNc = 13847 nodes. The names and 
codes of products are given in Table and their trade 
volumes, expressed in percent of the whole world trade 
volume, are given in Tablefor years 1998 and 2008. The 
main problem of construction of such a matrix is not its 
size, which is rather modest compared to those studied in 
uni, but the necessity to treat all countries on democratic 
grounds and at the same time to treat trade products on 
the basis of their trade volume. Indeed, the products can¬ 
not be considered on democratic grounds since their con¬ 


tributions to economy are linked with their trade volume. 
Thus, according to Table in year 2008 the trade volume 
of Petroleum and petroleum produets (code 33 in Table 
is by a factor 300 larger than those of Hides, skins and 
fur skins (undress.) (code 21 in Table f^. To incorporate 
these features in our mathematical analysis of multiprod¬ 
uct WTN we developed in this work the Google Person¬ 
alized Vector Method (GPVM) which allows to keep a 
democratic treatment of countries and at the same time 
to consider products proportionally to their trade volume. 
As a result we are able to perform analysis of the global 
multiproduct WTN keeping all interactions between all 
countries and all products. This is a new step in the WTN 
analysis since in our previous studies m it was possible 
to consider a trade between countries only in one product 
or only in all products summed together (all commodi¬ 
ties). The new finding of such global WTN analysis is an 
asymmetric ranking of products: some of them are more 
oriented to import and others are oriented to export while 
the ranking of products by the trade volume is always sym¬ 
metric after summation over all countries. This result with 
asymmetric ranking of products confirms the indications 
obtained on the basis ecological ranking [12] which also 
give an asymmetry of products in respect to import and 
export. Our approach also allows to analyze the sensitivity 
of trade network to price variations of a certain product. 

We think that the GPVM approach allows to perform a 
most advanced analysis of multiproduct world trade. The 
previous studies have been restricted to studies of statisti¬ 
cal characteristics of WTN links, patterns and their topol- 

PageRank algorithm to the WTN was discussed in [20] . 
the approach based on HITS algorithm was used in [21]. 
In comparison to the above studies, the approach devel¬ 
oped here for the multiproduct WTN has an advantage of 
analysis of ingoing and outgoing flows, related to PageR¬ 
ank and CheiRank, and of taking into account of multi¬ 
product aspects of the WTN. Even if the importance to 
multiproduct WTN analysis is clearly understood by re¬ 
searchers (see e.g. [22]) the Google matrix methods have 
not been efficiently used up to now. We also note that 
the matrix methods are extensively used for analysis of 
correlations of trade indexes (see e.g. umi]) but these 
matrices are Hermitian being qualitatively different from 
those appearing in the frame of Markov chains. Here we 
make the steps in multi-functional or multiproduct Google 
matrix analysis of the WTN extending the approach used 
in [TTj. 


2 Methods 


2.1 Google matrix construction for the WTN 

For a given year, we build Np money matrices Mf of the 
WTN from the COMTRADE database [T] (see [TT]). 

Mfc' = product p transfer (in USD) from country c' to c 

( 1 ) 
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Here the country indexes are c, c' = 1,..., A^c and a prod¬ 
uct index is p = 1,..., According to the COMTRADE 
database the number of UN registered countries is Nc = 
227 (in recent years) and the number of products is Np = 
10 and Np = 61 for 1 and 2 digits respectively from the 
Standard International Trade Classification (SITC) Rev. 
1. For convenience of future notation we also define the 
volume of imports and exports for a given country and 
product respectively as 

= = ( 2 ) 

d d 


The import and export volumes K* = 

shown for the world map of countries in Fig.f^ 
for year 2008. 

In order to compare later with PageRank and CheiRank 
probabilities we define volume trade ranks in the whole 
trade space of dimension N = Np x N^^ Thus the Impor- 
tRank (P) and ExportRank (P*) probabilities are given 
by the normalized import and export volumes 


Pi = VP/V, = 


( 3 ) 


where i = p + (c — l)A/p, i = 1 
volume is y = K,c' = 


,..., A^ and the total trade 

/-^p,C c Z-^p^C c 


The Google matrices G and G* are defined as N x N 
real matrices with non-negative elements: 


Gij = aSij^{l-a)viej , G*ij = a5'*ij-h(l-(r)u*ej , (4) 

where N = Np x Nc, a G {0^ 1] is the damping factor 
(0 < (a < 1), ej is the row vector of unit elements (e^ = 1), 
and Vi is a positive column vector called a personalization 
vector with Vi = 1 [ 7 ] • We note that the usual Google 
matrix is recovered for a personalization vector vi = ci/N 
In this work, following m, we fix a = 0.5. As discussed 
in laiiniiii] a variation of a in a range (0.5, 0.9) does not 
significantly affect the probability distributions of PageR¬ 
ank and CheiRank vectors. We specify the choice of the 
personalization vector a bit below. 

The matrices S and P* are built from money matrices 
M^cd as 


2.2 PageRank and CheiRank vectors from GPVM 


PageRank and CheiRank (P and P*) are defined as the 
right eigenvectors of G and G* matrices respectively at 
eigenvalue A = I: 

E = xd,Yl■ (6) 

j 3 


For the eigenstate at A = I we use the notation Pi = 
= 'ip'^i with the normalization YP — Yi^^^i — 1- 
For other eigenstates we use the normalization Yi \Pi?‘ — 
Yi \Pi\^ — 1- According to the Perron-Frobenius theorem 
the components of Pi, P*i are positive and give the prob¬ 
abilities to find a random surfer on a given node [7]. The 
PageRank K and CheiRank P* indexes are defined from 
the decreasing ordering of P and P* as P(P) > P(P + I) 
and P*(P) > P*(P* + I) with P, P* = I,..., P. 

If we want to compute the reduced PageRank and 
CheiRank probabilities of countries for all commodities 
(or equivalently all products) we trace over the product 
space getting = T,pPpc = Ep-P (P + (c - 1)-/Vp) and 
P* = EPpc = Y^pP* {p + {c-l)Np) with their corre- 
spending Pc and P* indexes. In a similar way we ob¬ 
tain the reduced PageRank and CheiRank probabilities 
for products tracing over all countries and getting 
Pp = Ec -P (^^ + (c - l)-y>) Ep Ppc and 
Pp = Yc {P Y Ppc with their correspond¬ 

ing product indexes Kp and P*. 

In summary we have Pp, P* = I,..., Np and Pc, P* = 
I,..., Pc- A similar definition of ranks from import and 
export trade volume can be done in a straightforward way 
via probabilities Pp, Pp , Pc, P^, Ppc^ Ppc correspond¬ 
ing indexes Kp, P*, Kc^ K^ , P, P*. 

To compute the PageRank and CheiRank probabilities 
from G and G* keeping democracy in countries and pro¬ 
portionality of products to their trade volume we use the 
GPVM approach with a personalized vector in Q. At the 
first iteration of Google matrix we take into account the 
relative product volume per country using the following 
personalization vectors for G and G*: 


VP 


Y*p 




■P' 


NcYp^Vc 




( 7 ) 


Si.., 




Kc'^pp/V:^ ifUc*"^0 
I/p if = 0 

Ml,Jp^p'/V!, ifUc^^O 
I/p if Vl, = 0 


( 5 ) 


where c,c' = I,... ,Nd, p,p' = 1,... ,Np]i =pP{c-l)Np] 
i' = p' P {c' — l)Np; and therefore i,i' = I,..., P. Note 
that the sum of each column of S and S'* are normalized 
to unity and hence the matrices G, G*, S, S* belong to the 
class of Google matrices and Markov chains. The eigen¬ 
values and eigenstates of G, G* are obtained by a direct 
numerical diagonalization using the standard numerical 
packages. 


using the definitions ([^ and the relation i = pP{c—l)Np. 
This personalized vector depends both on product and 
country indexes. In order to have the same value of person¬ 
alization vector in countries we can define the second it¬ 
eration vector proportional to the reduced PageRank and 
CheiRank vectors in products obtained from the GPVM 
Google matrix of the first iteration: 


v\i) 


Ne 





(8) 


In this way we keep democracy in countries but weighted 
products. This second iteration personalized vectors are 
used for the main part of computations and operations 
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with G and G*. This procedure with two iterations forms 
our GPVM approach. The difference between results ob¬ 
tained from the first and second iterations is not very 
large (see Figs.[^[^ but a detailed analysis of ranking of 
countries and products shows that the personalized vec¬ 
tor for the second iteration improves the results making 
them more stable and less fluctuating. In all Figures be¬ 
low (except Figs.j^l^ we show the results after the second 
iteration. 



Fig. 2. Dependence of probabilities of PageRank P{K), 
GheiRank ImportRank P{K) and ExportRank 

P*{K*) as a function of their indexes in logarithmic scale for 
WTN in 2008 with o = 0.5 at Nc = 227, Np = 1, N = 13847. 
Here the results for GPVM after and 2^^^ iterations are 
shown for PageRank (GheiRank) in red (blue) with dashed 
and solid curves respectively. ImportRank and ExportRank 
(trade volume) are shown by red and blue thin curves re¬ 
spectively. The ht exponents for PageRank and GheiRank are 
/3 = 0.61, 0.7 for the first iteration, /3 = 0.59, 0.65 for the second 
iteration, and 13 = 0.94,1.04 for ImportRank and ExportRank 
(for the range K G [10,2000]). 


The obtained results show the distribution of nodes 
on the PageRank-GheiRank plane (iC, i^*). In addition 
to two ranking indexes K, if* we use also 2DRank in¬ 
dex K 2 which combines the contribution of these indexes 
as described in [9]. The ranking list K 2 {i) is constructed 
by increasing K ^ K 3- I and increasing 2DRank in¬ 
dex K 2 {i) by one if a new entry is present in the list 
of first iC* < K entries of GheiRank, then the one unit 
step is done in R"* and K 2 is increased by one if the new 
entry is present in the list of first K < R"* entries of 
GheiRank. More formally, 2DRank K 2 {i) gives the order¬ 
ing of the sequence of sites, that appear inside the 
squares [1,1; K = if* = k; ...] when one runs progres¬ 
sively from k = 1 to N. Additionally, we analyze the dis¬ 
tribution of nodes for reduced indexes {Kp, (^c, )• 




Fig. 3. Probability distributions of PageRank and GheiRank 
for products Pp{Kp), Pp{Kp) (left panel) and countries 
Pc{Kc), iT(7^c) (right panel) in logarithmic scale for WTN 
from Fig. Here the results for the 1®^ and 2’^'^ GPVM itera¬ 
tions are shown by red (blue) curves for PageRank (GheiRank) 
with dashed and solid curves respectively. The probabilities 
from the trade volume ranking are shown by black curve (left) 
and dotted red and blue curves (right) for ImportRank and 
ExportRank respectively. 


We also characterize the localization properties of eigen¬ 
states of G, G* by the inverse participation ration (IPR) 
defined as ^ = (X^i This characteristic 

determines an effective number of nodes which contribute 
to a formation of a given eigenstate (see details in m)- 


2.3 Correlators of PageRank and GheiRank vectors 

Following previous works [Sl IOlITT] the correlator of PageR¬ 
ank and GheiRank vectors is defined as: 

N 

/t = Ar^P(i)P*(i) - 1. (9) 

i=l 

The typical values of k. are given in m for various net¬ 
works. 

For global PageRank and GheiRank the product-product 
correlator matrix is defined as: 

P(p^(c-l)Np)P^p'^{c-l)Np) 

Ec' - l)Np) P*{p' + (c" - l)Np)_ 

( 10 ) 



Then the correlator for a given product is obtained 
from (10) as: 

K,p — l^pp' 1 


where Sp^pf is the Kronecker delta. 

We also use the correlators obtained from the proba¬ 
bilities traced over products (Pc = E^p Ppc) arid over coun¬ 
tries {Pp = E^c Tpc) which are defined as 


Nc Np 

k{c) = iVe g PeP* - 1 , K{p) = iVp £ PpP; - 1 . (12) 

c=l p=l 

In the above equations (||)-@ the correlators are com¬ 
puted for PageRank and GheiRank probabilities. We can 
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also compute the same correlators using probabilities from 
the trade volume in ImportRank P and ExportRank P* 
defined by 

We discuss the values of these correlators in Section 4. 

3 Data description 

All data are obtained from the COMTRADE database 
[T]. We used products from COMTRADE SITC Rev. 1 
classification with number of products Np = 10 and 61. 
We choose SITC Rev.l since it covers the longest time 
interval. The main results are presented for Np = 61 with 
up to Nc = 227 countries. The names of products are 
given in Table their ImportRank index K and their 
fraction (in percent) of global trade volume in years 1998 
and 2008 are given in Table The data are collected 
and presented for the years 1962 - 2010. Our data and 
results are available at [25], the data for the matrices 
are available at COMTRADE [1] with the rules of their 
distribution policy. Following m we use for countries ISO 
3166-1 alpha-3 code available at Wikipedia. 

4 Results 


PageRank and CheiRank probabilities while there is no 
difference between ImportRank and ExportRank proba¬ 
bilities since they are equal after tracing over countries. 

After tracing over products we obtain probability dis¬ 
tributions Pc(Pc): Pci^c) countries shown in Fig.[^ 
We see that the probability of volume ranking varies ap¬ 
proximately by a factor 1000 while for PageRank and 
CheiRank such a factor is only approximate 10. Thus the 
democracy in countries induced by the Google matrix con¬ 
struction reduces significantly the variations of probabili¬ 
ties among countries and inequality between countries. 

Both panels of Fig. [^show relatively small variations 
between 1st and 2nd GPVM iterations confirming the sta¬ 
bility of this approach. In next sections we present the 
results only for 2nd GPVM iteration. This choice is con¬ 
firmed by consideration of ranking positions of various 
nodes of global matrices G, G* which show less fluctua¬ 
tions compared to the results of the 1st GPVM iteration. 

From the global ranking of countries and products we 
can select a given product and then determine local rank¬ 
ing of countries in a given product to see how strong is 
their trade for this product. The results for three selected 
products are discussed below for year 2008. For compar¬ 
ison we also present comparison with the export-import 
ranking from the trade volume. 


We apply the above methods to the described data sets of 
COMTRADE and present the obtained results below. 


4.1 PageRank and CheiRank probabilities 

The dependence of probabilities of PageRank P{K) and 
CheiRank P*{K'^) vectors on their indexes K, iC* are shown 
in Fig. I^for a selected year 2008. The results can be ap¬ 
proximately described by an algebraic dependence P oc 
1 /with the exponent values given in the caption. It is 
interesting to note that we find approximately the same 
P ^ 0.6 both for PageRank and CheiRank in contrast to 
the WWW, universities and Wikipedia networks where 
usually one finds ^ ^ 1 for PageRank and /3 ^ 0.6 for 
CheiRank lano]. We attribute this to an intrinsic property 
of WTN where the countries try to keep economy balance 
of their trade. The data show that the range of probability 
variation is reduced for the Google ranking compared to 
the volume ranking. This results from a democratic rank¬ 
ing of countries used in the Google matrix analysis that 
gives a reduction of richness dispersion between countries. 
The results also show that the variation of probabilities 
for 1st and 2nd GPVM results are not very large that 
demonstrates the convergence of this approach. 

After tracing probabilities over countries we obtain 
probability distributions Pp{Kp), Pp{Kp) over products 
shown in Fig. The variation range of probabilities is 
the same as for the case of volume ranking. This shows 
that the GPVM approach correctly treats products keep¬ 
ing their contributions proportional to their volume. The 
difference between 1st and 2nd iterations is rather small 
and is practically not visible on this plot. The impor¬ 
tant result well visible here is a visible difference between 


PageRank-CheiR^k 



ImportRank-ExportRank 



PageRank-CheiRank: all commodities 




. 1 1 , 

= ' * -i" : 


; i i j 

1 '" L. 

i 


-|" 4 -.- 



Fig. 4. Gountry positions on PageRank-GheiRank plane 
{Kc,K*) obtained by the GPVM analysis (left panels), 
ImportRank-ExportRank of trade volume (center panels), and 
for PageRank-GheiRank of all commodities (right panels, data 
from dH). Top panels show global scale {KcK* G [1,200]) 
and bottom panels show zoom on top ranks (i^c, N* ^ [1? 40])- 
Each country is shown by circle with its own flag (for a better 
visibility the circle center is slightly displaced from its integer 
position (Kc,Kc) along direction angle 7r/4). Data are shown 
for year 2008. 


4.2 Ranking of countries and products 

After tracing the probabilities P(7r), P*(i^*) over prod¬ 
ucts we obtain the distribution of world countries on the 
PageRank-CheiRank plane (Kc^K*) presented in Fig. 
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for a test year 2008. In the same figure we present the rank 
distributions obtained from Import Rank-Export Rank prob¬ 
abilities of trade volume and the results obtained in m 
for trade in all commodities. For the GPVM data we see 
the global features already discussed in HD: the countries 
are distributed in a vicinity of diagonal Kc = since 
each country aims to keep its trade balanced. The top 20 
list of top K 2 countries recover 15 of 19 countries of G20 
major world economies (EU is the number 20) thus ob¬ 
taining 79% of the whole list. This is close to the percent 
obtained in m for trade in all commodities. 

The global distributions of top countries with < 40, 

< 40 for the three ranking methods, shown in Fig. 
are similar on average. But some modifications introduced 
by the GPVM analysis are visible. Thus China (CHN) 
moves on 2nd position of CheiRank while it is in the 
1st position for trade volume ranking and CheiRank of 
all commodities. Also e.g. Saudi Arabia (SAU) and Rus¬ 
sia (RUS) move from the CheiRank positions i^*c = 21 
and K*c = 7 in all commodities E] to i^*c = 29 and 
= 6 in the GPVM ranking, respectively. Other ex¬ 
ample is a significant displacement of Nigeria (NGA). We 
explain such differences as the result of larger connectiv¬ 
ity required for getting high ranking in the multiproduct 
WTN. Indeed, China is more specialized in specific prod¬ 
ucts compared to USA (e.g. no petroleum production and 
export) that leads to its displacement in i^*c. We note 
that the ecological ranking gives also worse ranking po¬ 
sitions for China comparing to the trade volume ranking 
[ 12 ] . In a similar way the trade of Saudi Arabia is strongly 
dominated by petroleum and moreover its petroleum trade 
is strongly oriented on USA that makes its trade network 
concentrated on a few links while Russia is improving its 
position in due to significant trade links with EU and 
Asia. 

In global, the comparison of three ranks of countries 
shown in Fig. [^confirms that the GPVM analysis gives a 
reliable ranking of multiproduct WTN. Thus we now try 
to obtain new features of multiproduct WTN using the 
GPVM approach. 

The main new feature obtained within the GPVM ap¬ 
proach is shown in Fig. which gives the distribution of 
products on the Page Rank-CheiRank plane (Kp^Kp) af¬ 
ter tracing of global probabilities P{K), P*{K*) over all 
world countries. The data clearly show that the distribu¬ 
tion of products over this plane is asymmetric while the 
ranking of products from the trade volume produces the 
symmetric ranking of products located directly on diago¬ 
nal Kp = Kp. Thus the functions of products are asym¬ 
metric: some of them are more oriented to export (e.g. 
03 Fish and fish preparations, 05 Fruit and vegetables, 26 
Textile fibers, not manuf. etc., 28 Metalliferous ores and 
metal scrap, 84 Clothing); in last years (e.g. 2008) 34 Gas, 
natural and manufactured also takes well pronounced ex¬ 
port oriented feature characterized by location in the lower 
right triangle {Kf < Kp) of the square plane {Kp, K*). In 
contrast to that the products located in the upper left tri¬ 
angle {Kf > Kp) represent import oriented products (e.g. 
02 Dairy products and eggs, 04 Cereals and cereal prepara- 



Fig. 5. Two dimensional ranking of products on the 
PageRank-CheiRank plane {Kp,Kp). Each product is repre¬ 
sented by its specihc combination of color and symbol: color 
illustrates the hrst digit of COMTRADE SITC Rev. 1 code 
with the corresponding name shown in the legend on the right; 
symbols correspond to product names listed in Table with 
their code numbers; the order of names on the right panel of 
this Fig. is the same as in Table The trade volume ranking 
via ImportRank-ExportRank is shown by small symbols at the 
diagonal Kp = Kp, after tracing over countries this ranking is 
symmetric in products. Top left and right panels show years 
1963 and 1978, while bottom left and right panels show years 
1993 and 2008 respectively. 


tions, 64 Paper, paperboard and manuf., 65 Textile yarn, 
fabrics, etc., 86 Scientific & control instrum, photogr gds, 
clocks). 

It is interesting to note that the machinery products 
71, 72. 73 are located on leading import oriented positions 
in 1963, 1978, 1993 but they become more close to sym¬ 
metric positions in 2008. We attribute this to development 
of China that makes the trade in these products more sym¬ 
metric in import-export. It is interesting to note that in 
1993 the product 33 Petroleum and petroleum products 
loses its first trade volume position due to low petroleum 
prices but still it keeps the first CheiRank position show¬ 
ing its trade network importance for export. Each product 
moves on {Kp, Kf) with time. However, a part of the above 
points, we can say that the global distribution does not 
manifest drastic changes. Indeed, e.g. the green symbols 
of first digit 2 remain export oriented for the whole period 
1963 - 2008. We note that the established asymmetry of 
products orientation for the world trade is in agreement 
with the similar indications obtained on the basis of eco¬ 
logical ranking in [12]. However, the GPVM approach used 
here have more solid mathematical and statistical founda¬ 
tions with a reduced significance of fluctuations comparing 
to the ecological ranking. 

The comparison between the GPVM and trade vol¬ 
ume ranking methods provides interesting information. 
Thus in petroleum code 33 we have on top positions Rus¬ 
sia, Saudi Arabia, United Arab Emirates while from the 
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33 - Petroleum and petroleum products 72 - Electrical machinery, apparatus, and appliances 03 - Fish and fish preparations 



Fig. 6. Top panels show results of the GPVM data for coun¬ 
try positions on PageRank-GheiRank plane of local rank values 
K, K* ordered by {Kcp, K*p) for specific products with p = 33 
(left panel), p = 72 (center panel) and p = 03 (right panel). 
Bottom panels show the ImportRank-ExportRank planes re¬ 
spectively for comparison. Data are given for year 2008. Each 
country is shown by circle with its own flag as in Eig.|^ 

CheiRank order of this product we find Russia, USA, In¬ 
dia (see Fig. and Table [^. This marks the importance of 
the role of USA and India played in the WTN and in the 
redistribution of petroleum over nearby region countries, 
e.g. around India. Also Singapore is on a local petroleum 
position just behind India and just before Saudi Arabia, 
see Table This happens due to strong involvement of 
India and Singapore in the trade redistribution flows of 
petroleum while Saudi Arabia has rather restricted trade 
connections strongly oriented on USA and nearby coun¬ 
tries. 

For electrical machinery 72 there are less modifications 
in the top export or CheiRank positions (see Fig.|^ but we 
observe significant broadening of positions on PageRank- 
CheiRank plane comparing to ImportRank-ExportRank. 
Thus, Asian countries (China, Japan, S. Korea, Singa¬ 
pore) are located on the PageRank-CheiRank plane well 
below the diagonal K = showing a significant trade 
advantages of these countries in product 72 comparing to 
Western countries (USA, Germany, France, UK). 

Another product, shown in Fig. is 03 Fish and fish 
preparations. According to the trade volume export rank¬ 
ing the top three positions are attributed to China, Nor¬ 
way, Thailand. However, from CheiRank of product 03 we 
find another order with Thailand, USA, China. This result 
stresses again the broadness and robustness of the trade 
connections of Thailand and USA. As another example we 
note a significant improvement of Spain CheiRank posi¬ 
tion showing its strong commercial relations for product 
03. On the other side Russia has relatively good position 
in the trade volume export of 03 product but its CheiRank 
index becomes worse due to absence of broad commercial 
links for this product. 

The global top 20 positions of indexes K, i^*, 7^2, iC* 
are given in Table for year 2008. We note a signif¬ 
icant improvement of positions of Singapore and India 


in PageRank-CheiRank positions comparing to their po¬ 
sitions in the trade volume ranking. This reflects their 
strong commercial relations in the world trade. In the 
trade volume ranking the top positions are taken by 33 
petroleum and digit 7 of machinery products. This re¬ 
mains mainly true for PageRank-CheiRank positions but 
we see the spectacular improvement of positions of 84 
Clothing for China (TC* = 2) and 93 Special transact, 
for USA {K = 4) showing thus these two products have 
strong commercial exchange all over the world even if their 
trade volume is not so dominant. 



Fig. 7. Global plane of rank indexes (K,K*) for PageRank- 
GheiRank (left panel) and ImportRank-ExportRank (right 
panel) for N = 13847 nodes in year 2008. Each country and 
product pair is represented by a gray circle. Some countries 
are highlighted in colors: USA with black, South Korea with 
red, Ghina with green, Russia with red, France with yellow and 
Brazil with orange. 

We show the plane (7C, TC*) for the global world rank¬ 
ing in logarithmic scale in 2008 in Fig. The positions 
of trade nodes of certain selected countries are shown by 
color. We observe that the trade volume gives a higher 
concentration of nodes around diagonal comparing to the 
GPVM ranking. We attribute this to the symmetry of 
trade volume in products. 

In Fig. we show the distributions of top 200 ranks 
of the PageRank-GheiRank plane (zoom of left panel of 
Fig. Among the top 30 positions of TC* there are 8 
products of USA, 6 of China, 3 of Germany and other 
countries with less number of products. The top position 
at TC* = 1 corresponds to product 33 of Russia while Saudi 
Arabia is only at = 12 for this product. The lists of 
all N = 13847 network nodes with their 7C, K 2 ^ values 
are available at [25]. 

4.3 Time evolution of ranking 

The time evolution of indexes of products Kp, is shown 
in Fig. To obtain these data we trace PageRank and 
CheiRank probabilities over countries and show the time 
evolution of rank indexes of products Kp^K^ for top 15 
rank products of year 2010. The product 33 Petroleum 
and petroleum produets remains at the top CheiRank po¬ 
sition Kp = 1 for the whole period while in PageRank 
it shows significant variations from Kp = 1 to 4 being 
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Fig. 8. Top 200 global PageRank-CheiRank indexes (K, K*) 
distributions for year 2008. Each country (for different prod¬ 
ucts) is represented by its flag. 


at Kp = 4 at 1986 - 1999 when the petroleum had a low 
price. Products with first digit 7 have high ranks of Kp but 
especially strong variation is observed for Kp of 72 Elec¬ 
trical machinery moving from position 26 in 1962 to 4 in 
2010. Among other indexes with strong variations we note 
58 Plastic materials^ 84 Clothing^ 93 Special transact.^ 34 
Gas, natural and manufactured. 

The time evolution of products 33 and 72 on the global 
index plane {K^K"") is shown in Fig. 10 for 6 countries 
from Fig. Thus for product 72 we see a striking im¬ 
provement of iC* for China and S.Korea that is at the ori¬ 
gin of the global importance improvement of Kp in Fig.lol 
For the product 33 in Fig. pT| Russia improves significant^ 
its rank positions taking the top rank iC* = 1 (see also 
Table . 

The variation of global ranks K, if* with time is shown 
for 4 products and 10 countries in Fig. 11 For products 
72, 73 on a scale of 50 years we see a spectacular improve¬ 
ment of iC* for China, Japan, S.Korea. For the product 
33 we see strong improvement of iC* for Russia in last 15 
years. It is interesting to note that at the period 1986-1992 
of cheap petroleum 33 USA takes the top position if* = 1 
with a significant increase of its corresponding K value. 
We think that this is a result of political decision to make 
an economical pressure on USSR since such an increase of 
export of cheap price petroleum is not justified from the 
economical view point. For the product 33 we also note 
a notable improvement of iC* of India which is visible in 
CheiRank but not in ExportRank (see Table [^. We at¬ 
tribute this not to a large amount of trade volume but to 
a significant structural improvements of trade network of 
India in this product. We note that the strength and efff- 


K K* 

p p 



Fig. 9. Time evolution of PageRank Kp and CheiRank K* 
indexes for years 1962 to 2010 for certain products marked on 
the right panel side by their codes from Table Top panels 
show top 5 ranks of 2010, middle and bottom panels show 
ranks 6 to 10 and 11 to 15 for 2010 respectively. Colors of 
curves correspond to the colors of Fig. [^marking the first code 
digit. 



Fig. 10. Time evolution of ranking of two products 72 and 
33 for 6 countries of Fig. shown on the global PageRank- 
CheiRank plane (A, A*). Left and right panels show the cases 
of 72 Electrical machinery, apparatus and appliances and 33 
petroleum and petroleum products respectively. The evolution 
in time starts in 1962 (marked by cross) and ends in 2010 
(marked by square). 


ciency of trade network is also at the origin of significant 
improvement of PageRank and CheiRank positions of Sin¬ 
gapore comparing to the trade volume ranking. Thus the 
development of trade connections of certain countries sig¬ 
nificantly improves their Google rank positions. For the 
product 03 we note the improvement of A* positions of 
China and Argentina while Russia shows no improvements 
in this product trade for this time period. 
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Fig. 11. Time evolution of global ranking of PageRank and 
CheiRank indexes K* for selected 10 countries and 4 prod¬ 
ucts. Left and right panels show K and K* as a function of 
years for products: 03 Fish and fish preparations; 33 Petroleum 
and petroleum produets; 72, Eleetrieal maehinery, apparatus 
and applianees; and 73 Transport equipment (from top to bot¬ 
tom) . In all panels the ranks are shown in logarithmic scale for 
10 given countries: USA, South Korea, China, Russia, France, 
Brazil, Great Britain, Japan, Germany and Argentina marked 
by curve colors. 


4.4 Correlation properties of PageRank and CheiRank 


The properties of n correlator of PageRank and CheiRank 
vectors for various networks are reported in jSlITQ] . There 
are directed networks with small or even slightly negative 
values of e.g. Linux Kernel or Physical Review citation 
networks, or with /i: ^ 4 for Wikipedia networks and even 
larger values k ^ 116 for the Twitter network. 

The values of correlators defined by Eqs. are 

shown in Figs. 12|13 for a typical year 2008. For the global 
PageRank-CheiRank correlator we find k ^ 5.7 while 
for Import-Export probabilities the corresponding value is 
significantly larger with k. ^ 33.7. Thus the trade volume 
ranking with its symmetry in products gives an artificial 
increase of by a significant factor. A similar enhance¬ 
ment factor of Import-Export remains for correlators in 
products k{p) and countries k{c) from Eq. (12) while for 
PageRank-CheiRank we obtain a moderate correlator val¬ 
ues around unity (see Fig. 12). The PageRank-CheiRank 


correlator K,p 0 for specific products have relatively low 
values with Kp < 1 for practically all products with p < 45 
(we remind that here p counts the products in the order 
of their appearance in the Table it is different from 
COMTRADE code number). 

The correlation matrix of products hippf (10) is shown 
in Fig. This matrix is asymmetric and demonstrates 
the existence of relatively high correlations between prod¬ 
ucts 73 Transport equipment^ 65 Textile yarn, fabries, made 



Fig. 12. PageRank-CheiRank correlators np (11) from the 


GPVM are shown as a function of the product index p with 
the corresponding symbol from Fig. PageRank-CheiRank 
and Import Rank-Export Rank correlators are shown by solid 
and dashed lines respectively, where the global correlator n 
it is shown in black, the correlator for countries k ,{ c ) 
is shown by red lines, the correlator for products K,{p) is 
shown by blue lines. Here product number p is counted in order 
of appearance in Table^ The data are given for year 2008 with 
iVp = 61, iVc = 227, N = 13847. 



012345 6789 
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Fig. 13. Product PageRank-CheiRank correlation matrix 
tvp^pf (10) for year 2008 with correlator values shown by color. 
The code indexes p and p of all Np = 61 products are shown 
on X and y axes by their corresponding first digit (see Table [^. 


up artieles, ete. and 83 Travel goods, handbags and similar 
artieles that all are related with transportation of prod¬ 
ucts. 


4.5 Spectrum and eigenstates of WTN Google matrix 

Above we analyzed the properties of eigenstates of G and 
G* at the largest eigenvalue A = 1. However, in total there 
are N eigenvalues and eigenstates. The results obtained 
for the Wikipedia network [26] demonstrated that eigen¬ 
states with large modulus of A correspond to certain spe¬ 
cific communities of the network. Thus it is interesting to 
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study the spectral properties of G for the multiproduct 
WTN. The spectra of G and G* are shown in Fig, for 
year 2008. It is interesting to note that for G the spectrum 
shows some similarities with those of Wikipedia (see Fig.l 
in [26]). At (a = 1 there are 12 and 7 degenerate eigenval¬ 
ues A = 1 for G and G* respectively. Thus the spectral 
gap appears only for a < 1. The dependence of IPR ^ of 
eigenstates of G on ReX is shown in Fig. The results 
show that ^ ^ so that the eigenstates are well localized 
on a certain group on nodes. 



Fig. 14. Spectrum of Google matrices G (left panel) and G* 
(right panel) represented in the complex plane of A. The data 
are for year 2008 with 0 = 1, and N = 13847, Nc = 227, 
Np = 61. Four eigenvalues marked by colored circles are used 
for illustration of eigenstates in Figs. Eau 


The eigenstates can be ordered by their decreas¬ 
ing amplitude \ giving the eigenstate index Ki with the 
largest amplitude at = 1. The examples of four eigen¬ 
states are shown in Fig. We see that the amplitude is 
mainly localized on a few top nodes in agreement of small 
values of ^ ^ 4 shown in Fig. The top ten amplitudes 
of these four eigenstates are shown in Table with cor¬ 
responding names of countries and products. We see that 
for a given eigenstate these top ten nodes correspond to 
one product clearly indicating strong links of trade be¬ 
tween certain countries. Thus for 06 Sugar we see strong 
link between geographically close Mali and Guinea with 
further links to USA, Germany etc. In a similar way for 56 
Fertilizers there is a groups of Latin American countries 
Brazil, Bolivia, Paraguay linked to Argentina, Uruguay 
etc. We see a similar situation for products 57 and 52. 
These results confirm the observation established in [26] 
for Wikipedia that the eigenstates with large modulus of A 
select interesting specific network communities. We think 
that it would be interesting to investigate the properties 
of eigenstates in further studies. 


4.6 Sensitivity to price variations 

Above we established the global mathematical structure 
of multiproduct WTN and presented results on its ranking 
and spectral properties. Such ranking properties bring new 
interesting and important information about the WTN. 
However, from the view point of economy it is more im¬ 
portant to analyze the effects of crisis contamination and 
price variations. Such an analysis represents a complex 


80 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

W 

60- ‘ ■ 



Re(A) 


Fig. 15. Inverse participation ratio (IPR) ^ of all eigenstates 
of G as a function of the real part of the corresponding eigen¬ 
value A from the spectrum of Fig. [7^ T he eigenvalues marked 
by color circles are those from Fig. |14| 



Fig. 16. Eigenstate amplitudes ordered by its own de¬ 
creasing amplitude order with index Ki for 4 different eigen¬ 
values of Fig. 14 (states are normalized as |'0i| = 1). Top 


panel shows two example of real eigenvalues with A = 0.9548 
and A = 0.9345 while bottom panel shows two eigenval¬ 
ues with large imaginary part with A = 0.452 -h iO.775 and 
A = 0.424-hiO.467. Node names (country, product) for top ten 
largest amplitudes of these eigenvectors are shown in Table ^ 


task to which we hope to return in our further investiga¬ 
tions. However, the knowledge of the global WTN struc¬ 
ture is an essential building block of this task and we think 
that the presented results demonstrate that this block is 
available now. 

Using the knowledge of WTN structure, we illustrate 
here that it allows to obtain nontrivial results on sensitiv¬ 
ity to price variations for certain products. We consider as 
an example year 2008 and assume that the price of prod¬ 
uct 33 Petroleum and petroleum produets is increased by 
a relative fraction S going from its unit value 1 to 1 + 
(or 6 = (^ 33 ). Then we compute the derivatives of proba- 
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Fig. 17. Derivatives D — dP/dSss and D* = dP"^/dSss for a 
price variation S 33 of 33 Petroleum and petroleum produets” for 
year 2008. Top left and right panels show the cases of negative 
and positive D and D* respectively, with some products and 
countries labeled by their 2 digit code. Bottom panels show 
the positive and negative cases of the logarithmic derivatives 
Di = D/P and D/ = D""/P* for countries and products with 
K 2 < 50, where the flags and 2 digit codes for countries and 
products are shown (in right panels only product 33 is present). 
Codes are described in Table [T] 


bilities of PageRank D = dPjdd = AP /6 and CheiRank 
1)* = dP^'/dd = AP'^/S. The computation is done for 
values of (5 = 0.01,0.03,0.05 ensuring that the result is 
not sensitive to a specific S value. We also compute the 
logarithmic derivatives Di = dlnP/dS, D/ = dlnP^/dS 
which give us a relative changes of P, P*. 
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Fig. 18. Same as in top panels of Fig. but using probabil¬ 
ities from the trade volume §. 


and For (P,P*) the nodes are distributed in 

two sectors with P > 0, P* > 0 and D < <0. The 

largest values with P > 0, P* > 0 correspond to nodes of 
countries of product 33 which are rich in petroleum (e.g. 
Russia, Saudi Arabia, Nigeria) or those which have strong 
trade transfer of petroleum to other countries (Singapore, 
India, China etc). It is rather natural that with the growth 
of petroleum prices the rank probabilities P, P* of these 
countries grow. A more unexpected effect is observed in 
the sector P < 0, P* < 0. Here we see that an increase 
of petroleum price leads to a decrease of probabilities of 
nodes of countries Germany, France, China, Japan trading 
in machinery products 71, 72, 73. 

For comparison we also compute the derivatives P, P*, 
Di^D/ from the probabilities (§ defined by the trade 
volume of Import-Export insteacTof PageRank-CheiRank. 
The results are shown in Fig. for petroleum price vari¬ 
ation to be compared with Fig. The distribution of 
P, P* is rather different from those values obtained with 
PageRank-CheiRank probabilities. This is related to the 
fact that PageRank and CheiRank take into account the 
global network structure while the trade volume gives only 
local relations in trade links between countries. The dif¬ 
ference between these two methods becomes even more 
striking for logarithmic derivatives P/, D/. Indeed, for the 
trade volume ranking the variation of probabilities P*, P 
due to price variation of a given product can be com¬ 
puted analytically taking into account the trade volume 
change with 6 p. The computations give Dcp = (1 —/p)Pcp, 
D\p = (1 - fp)P^ 


cp for a derivative of of probability of 

and 


product p and country c over the price of product Sp 


P, 


cp' = -fpPcp', Dcp' = -fpPcp', D*cp' = -fpP*cp' (if 
p' 7^ p); where fp is a fraction of product p in the world 
trade. From these expressions we see that the logarithmic 
derivatives are independent of country and product. In¬ 
deed, for the case of Fig. we obtain analytically and by 
direct numerical computations that Di = P* = —0.2022 
(for all countries if p' ^ p = 33) and Di = D/ = 0.7916 
(for all countries if p' = p = 33). Due to simplicity of this 
case we do not show it in Fig. 

The results for price variation of 34 Gas, natural and 
manufactured are presented in Fig, showing derivatives 
of PageRank and CheiRank probabilities over J 34 . We see 
that for absolute derivatives P, P* the mostly affected 
are now nodes of gas producing countries for the sector 
P, P* > 0, while for the sector P, P* < 0 the mostly 
affected are countries linked to petroleum production or 
trade, plus USA with products 71,72,73. For the sector 
of logarithmic derivatives Di , D/ < 0 among top K 2 and 
K, nodes we find nodes of countries of product 33 and 
also 93. 

Thus the analysis of derivatives provides an interest¬ 
ing new information of sensitivity of world trade to price 
variations. 


The results for the price variation J 33 of 33 Petroleum 
and petroleum products are shown in Fig. The deriva¬ 
tives for all WTN nodes are shown on the planes (P, P*) 


4.7 World map of CheiRank-PageRank trade balance 

On the basis of the obtained WTN Google matrix we can 
now analyze the trade balance in various products between 
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Di D, 

Fig. 19. Derivative of P and P* (D and P* respectively) for 
a price variation of 34 Gas, natural and manufactured ioi 2008. 
Top left and right panels show the cases of negative and posi¬ 
tive sectors of D and P* respectively, with some products and 
countries labeled by their 2 digit code and names (in top right 
panel all points correspond to product 34). Bottom panels show 
the cases of the logarithmic derivatives Di and D* for countries 
and products with K 2 < 50 (bottom left panel) and K, < 
25 (bottom right panel); flags and 2 digit codes for countries 
and products are shown. In bottom right panel (P, K* < 25) 
we do not show the case of Sudan {73 Transport equipment) 
which has values of {Di,Di) = (2 x 10“^, 1.75 x 10“^). Codes 
are described in Table [T] 


the world countries. Usually economists consider the ex¬ 
port and import of a given country as it is shown in Fig.[2 
Then the trade balance of a given country c can be defined 
making summation over all products: 

Be = = (p;-Pc)/(p;+^c). 

p p 

(13) 

In economy, Pc^Pq defined via the probabilities of 
trade volume Pcp^ pk. In our approach, we define 

Pep, Pep as PageRank and CmeiRank probabilities. In con¬ 
trast to the trade volume our approach takes into account 
the multiple network links between nodes. 

The comparison of the world trade balance obtained 
by these two methods is shown in Fig.[^ We see that the 
leadership of China becomes very well visible in CheiRank- 
PageRank balance map while it is much less pronounced 
in the trade volume balance. The Google matrix analysis 
also highlights the dis-balance of trade network of Nige¬ 
ria (strongly oriented on petroleum export and machinery 
import) and Sudan. It is interesting to note that the posi¬ 
tive CheiRank-PageRank balance is mainly located in the 
countries of BRIGS (Brazil, Russia, India, China, South 
Africa). In contrast to that, the usual trade volume bal¬ 



Fig. 20. World map of probabilities balance Be = {Pf — 
Pc) / (P* +Pc) determined for each of Nc = 227 countries in year 
2008. Top panel: probabilities Pf, Pc are given by CheiRank 
and PageRank vectors; bottom panel: probabilities are com¬ 
puted from the trade volume of Export-Import Names of 
countries can be find at |4]. 


ance highlights Western Sahara and Afganistan at large 
positive and negative trade balance in 2008. 



Fig. 21. Derivative of probabilities balance dBc/dSss over 
petroleum price S33 for year 2008. Top panel: balance of coun¬ 
tries Be is determined from CheiRank and PageRank vectors 
as in the top panel of Fig. bottom panel: Be values are 
computed from the trade volume as in the bottom panel of 
Fig. 20 Names of countries can be find at |4]. 
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We can also determine the sensitivity of trade balance 
to price variation of a certain product p computing the 
balance derivative dBc/d 6 p. The world map sensitivity in 
respect to price of petroleum p = 33 is shown in Fig. for 
the above two methods of definition of probabilities Pc, Pq 
in (13). For the CheiRank-PageRank balance we see that 
the derivative dBdddss is positive for countries produc¬ 
ing petroleum (Russia, Saudi Arabia, Venezuela) while 
the highest negative derivative appears for China which 
economy is happened to be very sensitive to petroleum 
price. The results from the trade volume computation of 
dBc/d 6 p^ shown in Fig. 21, give rather different distribu¬ 


tion of derivatives over countries with maximum for Su¬ 
dan and minimum for the Republic of Nauru (this coun¬ 
try has very small area and is not visible in the bottom 
panel of Fig. [21] ), while for China the balance looks to be 
not very sensitive to 6 ss (in contrast to the CheiRank- 
PageRank method). This happens due to absence of links 
between nodes in the trade volume computations while the 
CheiRank-PageRank approach takes links into account 
and recover hidden trade relations between products and 
countries. 



Fig. 22. Derivative of partial probability balance of product 
p dehned as dBcpjdSss over petroleum price S 33 for year 2008; 
here Bcp = {Pep — Pep)/{Pc + Pc) and p — 72 {72 Electrical 
machinery ... from Table [^; the product balance of countries 
Bcp is determined from CheiRank and PageRank vectors (top 
panel) and from the trade volume of Export-Import it (bot¬ 
tom panel). Names of countries can be hnd at [4]. 


This absence of links in the trade volume approach 
becomes also evident if we consider the derivative of the 
partial trade balance for a given product p defined as 

= {p*^-p,^)/j2ip:p+Pcp) = {p:^-Pcp)/{p:+Pc), 

p 

(14) 


SO that the global country balance is Be = ^p Bcp. Then 
the sensitivity of partial balance of a given product p in 
respect to a price variation of a product p' is given by 
the derivative dBcp/ddp'. The sensitivity for balance of 
product p = 12 {72 Electrical machinery ...) in respect 
to petroleum p' = 33 price variation ^33 is shown for the 
CheiRank-PageRank balance in Fig. [^ (top panel) in¬ 
dicating sensitivity of trade balance of product p = 72 
at the petroleum p' = 33 price variation. We see that 
China has a negative derivative for this partial balance. 
In contrast, the computations based on the trade volume 
bottom panel) give a rather different distribu- 
erNatives dBcp/d 6 p> over countries. In the trade 
volume approach the derivative dBcp/d 6 p> appears due to 
the renormalization of total trade volume and nonlinearity 
coming from the ratio of probabilities. We argue that the 
CheiRank-PageRank approach treats the trade relations 
between products and countries on a significantly more 
advanced level taking into account all the complexity of 
links in the multiproduct world trade. 

Using the CheiRank-PageRank approach we determine 
the sensitivity of partial balance of all 61 products in re¬ 
spect to petroleum price variation ^33 for China, Russia 
and USA, as shown in Fig. (top panel). We see that 
the diagonal derivative dBcsT/d^^^ positive for Russia 
but is negative for China and USA. Even if USA produce 
petroleum its sensitivity is negative due to a significant 
import of petroleum to USA. For non-diagonal derivatives 
over 6 s 3 we find positive sensitivity of Russia and USA 
for products p = 71, 72, 73 while for China it is negative. 
Other product partial balances sensitive to petroleum are 
e.g. 84 Clothing for China for which expensive petroleum 
gives an increase of transportation costs; negative deriva¬ 
tive of balance in metal products p = 67, 68 for Russia 
due to fuel price increase; positive derivative for 93 Spe¬ 
cial transact. ... of USA. 


(Fig. [^ 
tion of c 


The sensitivity of country balance Be to price vari¬ 
ation Spf for all products is shown in Fig. [^ for China 
(middle panel) and USA (bottom panel). We find that 
the balance of China is very sensitive to p' = 33, 84 and 
indeed, these products play an important role in its econ¬ 
omy with negative and positive derivatives respectively. 
For USA the trade balance is also very sensitive to these 
two products p' = 33, 84 but the derivative is negative 
in both cases. We also present the derivative of balance 
without diagonal term {d{Bc — Bcp')/d 6 p') for China and 
USA. This quantity shows that for USA all other products 
give a positive derivative for p' = 33 but the contribution 
of petroleum import gives the global negative derivative 
of the total USA balance. In a similar way for China for 
p' = 84 all products, except the diagonal one p' = 84, give 
a negative sensitivity for balance but the diagonal con¬ 
tribution of p' = 84 gives the final positive derivative of 
China total balance in respect to 6 sa. 


The CheiRank-PageRank approach allows to deter¬ 
mine cross-product sensitivity of partial trade balance com¬ 
puting the derivative dBcp/ddpi shown in Fig. 24 for China 
and USA. The derivatives are very different for two coun¬ 
tries showing a structural difference of their economies. 
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Fig. 23. Top panel: derivative dBcp/dSss of partial probability 
balance Bcp of product p over petroleum price S33 for year 2008 
and countries: China (black circles), Russia (blue squares) and 
USA (red diamonds); inset panel shows the products of digit 
3 including the diagonal term p = 33 being out of scale in 
the main panel; here Bcp = {Pep — Pep)/{Pc P Pc) |T^ . Center 
(China) and bottom (USA) panels show derivative dBc/d6p/ of 
country total probability balance Be over price 5^/ of product 
p' for year 2008; derivatives of balance without diagonal term 
{dBc/dSpf — dBcp'/ddpi) are represented by open circles and 
open diamonds for China and USA respectively. The product 
balance of countries Bcp and Be are determined from CheiRank 
and PageRank vectors. The vertical dotted lines mark the first 
digit of product index p or p from Table [l] 


Thus for China the cross-derivative (at p ^ p') are mainly 
negative (except a few lines around p = 33) but the diag¬ 
onal terms dBcp/d6p are mainly positive. In contrast, for 
USA the situation is almost the opposite. We attribute 
this to the leading role of China in export and the lead¬ 
ing role of USA in import. However, a detailed analysis of 
these cross-products derivatives and correlations require 
further more detailed analysis. We think that the pre¬ 
sented cross-product sensitivity plays and important role 
in the multiproduct trade network that are highlighted by 
the Google matrix analysis developed here. This analysis 
allows to determine efficiently the sensitivity of multiprod¬ 
uct trade in respect to price variations of various products. 


China 



0.01 

0.005 

0 


-0.005 


- 0.01 


0.01 

0.005 

0 


-0.005 


- 0.01 


Fig. 24. China (top) and USA (bottom) examples of deriva¬ 
tive dBcp/dSpf of partial probability balance Bcp of product 
p over price Spf of product p' for year 2008 . Diagonal terms, 
given by dBcp/ddp ys. p — p ^ are shown on the top panels of 
each example. Products p and p are shown in x-axis and y-axis 
respectively (indexed as in Table [R, while dBcp/d6p/ is repre¬ 
sented by colors with a threshold value given by —0.01 and 
0.01 for negative and positive values respectively, also shown 
in red dashed lines on top panels with diagonal terms. Dotted 
lines mark the first digit of Table Here Bcp are defined by 
CheiRank and PageRank probabilities. 


5 Discussion 

In this work we have developed the Google matrix analysis 
of the multiproduct world trade network. Our approach 
allows to treat all world countries on equal democratic 
grounds independently of their richness keeping the con¬ 
tributions of trade products proportional to their fractions 
in the world trade. As a result of this approach we have 
obtained a reliable ranking of world countries and prod¬ 
ucts for years 1962 - 2010. The Google analysis captures 
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the years with crises and also shows that after averaging 
over all world countries some products are export oriented 
while others are import oriented. This feature is absent in 
the usual Import-Export analysis based on trade volume 
which gives a symmetric orientation of products after such 
an averaging. 


The WTN matrix analysis determines the trade bal¬ 
ance for each country not only in trade volume but also in 
CheiRank-PageRank probabilities which take into account 
multiple trade links between countries which are absent in 
the usual Export-Import considerations. The CheiRank- 
PageRank balance highlights in a clear manner the lead¬ 
ing WTN role of new rising economies of China and other 
BRICS countries. This analysis also allows to determine 
the sensitivity of trade network to price variations of var¬ 
ious products that opens new possibilities for analysis of 
cross-product price influence via network links absent in 
the standard Export-Import analysis. 


We think that this work makes only first steps in the 
development of WTN matrix analysis of multiproduct world 
trade. Indeed, the global properties of the Google matrix 
of multiproduct WTN should be studied in more detail 
since the statistical properties of matrix elements of G, 
shown in Fig. for year 2008, are still not well under¬ 
stood (e.g. visible patterns present in the coarse-grained 
representation of G in Fig. 


23). 


Even if the UN COMTRADE database contains a lot 
of information there are still open questions if all essential 
economic aspects are completely captured in this database. 
Indeed, the COMTRADE data for trade exchange are di¬ 
agonal in products since there are no interactions (trade) 
between products. However, this feature may be a weak 
point of collected data since in a real economy there is 
a transformation of some products into some other prod¬ 
ucts (e.g. metal and plastic are transferred to cars and 
machinery). It is possible that additional data should be 
collected to take into account the existing interactions be¬ 
tween products. There are also some other aspects of ser¬ 
vices and various other activities which are not present 
in the COMTRADE database and which can affect the 
world economy. 


One of the important missing element of COMTRADE 
are financial flows between co untri es. Indeed, the product 
93 Special trans. ... (see Tables [m> partially takes into ac¬ 
count the financial flows but it is clear that the interbank 
flows are not completely reported in the database. In fact 
the Wold Bank Web (WBW) really exists (e.g. a private 
person can transfer money from his bank account to an¬ 
other person account using SWIFT code) but the flows on 
the WBW remain completely hidden and not available for 
scientific analysis. The size on interbank networks are rel¬ 
atively small (e.g. the whole Federal Reserve of USA has 
only N ^ 6600 bank nodes m and there are only about 
N ^ 2000 bank nodes in Germany [28]). Thus the WBW 
size of the whole world is about a few tens of thousands 
of nodes and the Google matrix analysis should be well 
adapted for WBW. We consider that there are many sim¬ 
ilarities between the multiproduct WTN and the WBW, 
where financial transfers are performed with various finan¬ 


cial products so that the above WTN analysis should be 
well suited for the WBW. The network approach to the 
WBW flows is now at the initial development stage (see 
e.g. [271[28l|29]) but hopefully the security aspects will be 
handled in an efficient manner opening possibilities for the 
Google matrix analysis of the WBW. The joint analysis of 
trade and financial flows between world countries would 
allow to reach a scientific understanding of peculiarities 
of such network flows and to control in an efficient way 
financial and petroleum crises. 

The developed Google matrix analysis of multiprod¬ 
uct world trade allows to establish hidden dependencies 
between various products and countries and opens new 
prospects for further studies of this interesting complex 
system of world importance. 
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Fig. 25. Google matrix Gkk' representation for 2008 with 
a = 0.5 ordered by PageRank index iLvalue (where K — K' — 
1 is on top left corner). Top panel shows the whole Google 
matrix [N — Nc'X Np = 227 x 61 = 13847) with coarse-graining 
of iV X A" elements down to 200 x 200 shown cells. Genter 
panel represents the top corner of the full Google matrix with 
K,K' < 200. Bottom panel shows the coarse-grained Google 
matrix for countries for the top 100 countries {Kc,Kc < 100). 
Golor changes from black at minimal matrix element to white 
at maximal element, a = 0.5. 
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Table 1. Codes and names of the 61 products from COMTRADE Standard International Trade Classihcation (SITC) Rev. 1. 


code 

name 

code 

name 

00 

Live animals 

54 

Medicinal and pharmaceutical products 

01 

Meat and meat preparations 

55 

Perfume materials, toilet & cleansing preptions 

02 

Dairy products and eggs 

56 

Eertilizers, manufactured 

03 

Eish and fish preparations 

57 

Explosives and pyrotechnic products 

04 

Cereals and cereal preparations 

58 

Plastic materials, etc. 

05 

Emit and vegetables 

59 

Chemical materials and products, nes 

06 

Sugar, sugar preparations and honey 

61 

Leather, Ithr. Manufs., nes Sz dressed fur skins 

07 

Coffee, tea, cocoa, spices & manufacs. Thereof 

62 

Rubber manufactures, nes 

08 

Eeed. Stuff for animals excl. Unmilled cereals 

63 

Wood and cork manufactures excluding furniture 

09 

Miscellaneous food preparations 

64 

Paper, paperboard and manufactures thereof 

11 

Beverages 

65 

Textile yarn, fabrics, made up articles, etc. 

12 

Tobacco and tobacco manufactures 

66 

Non metallic mineral manufactures, nes 

21 

Hides, skins and fur skins, undressed 

67 

Iron and steel 

22 

Oil seeds, oil nuts and oil kernels 

68 

Non ferrous metals 

23 

Crude rubber including synthetic and reclaimed 

69 

Manufactures of metal, nes 

24 

Wood, lumber and cork 

71 

Machinery, other than electric 

25 

Pulp and paper 

72 

Electrical machinery, apparatus and appliances 

26 

Textile hbres, not manufactured, and waste 

73 

Transport equipment 

27 

Crude fertilizers and crude minerals, nes 

81 

Sanitary, plumbing, heating and lighting hxt. 

28 

Metalliferous ores and metal scrap 

82 

Eurniture 

29 

Crude animal and vegetable materials, nes 

83 

Travel goods, handbags and similar articles 

32 

Coal, coke and briquettes 

84 

Clothing 

33 

Petroleum and petroleum products 

85 

Eootwear 

34 

Gas, natural and manufactured 

86 

Scientif Sz control instrum, photogr gds, clocks 

35 

Electric energy 

89 

Miscellaneous manufactured articles, nes 

41 

Animal oils and fats 

91 

Postal packages not class. According to kind 

42 

Eixed vegetable oils and fats 

93 

Special transact. Not class. According to kind 

43 

Animal and vegetable oils and fats, processed 

94 

Animals, nes, inch Zoo animals, dogs and cats 

51 

Chemical elements and compounds 

95 

Eirearms of war and ammunition therefor 

52 

53 

Crude chemicals from coal, petroleum and gas 
Dyeing, tanning and colouring materials 

96 

Coin, other than gold coin, not legal tender 
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Table 2. Columns represent data: codes of 61 products of COMTRADE SITC Rev.l, ImportRank and ExportRank K = K* 
in year 2008, product fraction in global trade volume in 2008, K = K* in 1998, product fraction in 1998. 


code 

K(08) 

% vol(08) 

is:(98) 

% vol (98) 

code 

K(08) 

% vol(08) 

iC(98) 

% vol (98) 

00 

53 

0.10 

51 

0.17 

54 

9 

2.89 

16 

1.88 

01 

27 

0.69 

26 

0.83 

55 

25 

0.76 

28 

0.79 

02 

34 

0.44 

34 

0.56 

56 

30 

0.55 

43 

0.36 

03 

28 

0.63 

22 

0.99 

57 

58 

0.03 

57 

0.04 

04 

21 

1.07 

19 

1.13 

58 

15 

1.95 

13 

2.07 

05 

19 

1.16 

18 

1.50 

59 

22 

1.04 

20 

1.13 

06 

49 

0.23 

44 

0.36 

61 

51 

0.19 

42 

0.37 

07 

33 

0.47 

29 

0.73 

62 

26 

0.73 

24 

0.85 

08 

38 

0.39 

36 

0.45 

63 

35 

0.43 

32 

0.61 

09 

40 

0.34 

41 

0.39 

64 

20 

1.14 

17 

1.79 

11 

31 

0.54 

31 

0.65 

65 

18 

1.40 

11 

2.46 

12 

47 

0.24 

35 

0.49 

66 

17 

1.71 

14 

2.01 

21 

56 

0.05 

53 

0.11 

67 

7 

3.63 

10 

2.74 

22 

37 

0.39 

48 

0.32 

68 

11 

2.27 

15 

1.95 

23 

44 

0.26 

50 

0.22 

69 

13 

2.04 

12 

2.12 

24 

39 

0.35 

30 

0.65 

71 

2 

11.82 

1 

15.03 

25 

43 

0.29 

45 

0.34 

72 

3 

10.42 

3 

12.26 

26 

50 

0.22 

37 

0.45 

73 

4 

10.06 

2 

12.38 

27 

41 

0.33 

47 

0.33 

81 

42 

0.31 

46 

0.34 

28 

16 

1.92 

25 

0.84 

82 

23 

0.93 

21 

1.03 

29 

48 

0.24 

40 

0.39 

83 

45 

0.26 

49 

0.26 

32 

24 

0.82 

39 

0.42 

84 

10 

2.42 

6 

3.44 

33 

1 

14.88 

4 

5.02 

85 

29 

0.59 

27 

0.79 

34 

14 

2.04 

23 

0.99 

86 

12 

2.25 

8 

2.95 

35 

46 

0.26 

52 

0.17 

89 

6 

3.72 

5 

4.54 

41 

57 

0.03 

58 

0.04 

91 

61 

0.00 

61 

0.00 

42 

32 

0.49 

38 

0.44 

93 

5 

3.92 

9 

2.92 

43 

54 

0.08 

55 

0.08 

94 

59 

0.01 

59 

0.01 

51 

8 

3.01 

7 

3.07 

95 

55 

0.08 

54 

0.11 

52 

52 

0.11 

56 

0.05 

96 

60 

0.00 

60 

0.00 

53 

36 

0.39 

33 

0.60 







Table 3. Top 20 ranks for global PageRank K, CheiRank K*, 2dRank K 2 , ImportRank K and ExportRank K* for given 
country and product code for year 2008. 


# 

K 

country Sz code 

K- 

country Sz code 

K 2 

country Sz code 

K 

country & code 

k* 

country Sz code 

1 

USA 

33 

Russia 

33 

Germany 

73 

USA 

33 

China 

72 

2 

USA 

73 

China 

84 

USA 

73 

USA 

71 

Russia 

33 

3 

USA 

71 

Germany 

73 

USA 

33 

USA 

72 

China 

71 

4 

USA 

93 

Japan 

73 

USA 

71 

USA 

73 

Germany 

73 

5 

Germany 

73 

USA 

73 

India 

33 

Japan 

33 

Germany 

71 

6 

USA 

72 

China 

72 

Singapore 

33 

China 

72 

Saudi Arabia 

33 

7 

Erance 

73 

USA 

33 

Germany 

71 

China 

33 

USA 

71 

8 

Germany 

71 

India 

33 

USA 

72 

Germany 

71 

Japan 

73 

9 

Singapore 

33 

USA 

71 

Erance 

73 

Germany 

73 

USA 

73 

10 

India 

33 

China 

71 

Netherlands 

33 

Netherlands 

33 

Japan 

71 

11 

China 

33 

Singapore 

33 

USA 

93 

Germany 

72 

USA 

72 

12 

Netherlands 

33 

Saudi Arabia 

33 

Nigeria 

33 

China 

71 

China 

89 

13 

Erance 

33 

Germany 

71 

Germany 

72 

USA 

89 

Germany 

72 

14 

UK 

71 

USA 

72 

China 

72 

Italy 

33 

China 

84 

15 

UK 

73 

Erance 

73 

China 

71 

Germany 

33 

Japan 

72 

16 

Germany 

72 

Thailand 

3 

UK 

33 

South Korea 

33 

South Korea 

72 

17 

USA 

89 

Kazakhstan 

33 

Germany 

93 

Erance 

73 

Erance 

73 

18 

South Korea 

33 

U. Arab Emir. 

33 

China 

33 

China 

28 

Italy 

71 

19 

Erance 

71 

USA 

28 

South Korea 

33 

Germany 

93 

U. Arab Emir. 

33 

20 

Sudan 

73 

Netherlands 

33 

Australia 

33 

India 

33 

Germany 

93 
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Table 4. Top 10 values of 4 different eigenvectors from Fig. 16 The corresponding eigenvalues form left to right are A = 0.9548, 


A = 0.9345, A = 0.452+ z0.775 and A = 0.424+ z0.467. There is only one product in each of these top 10 list nodes which are: 57 
Explosives and pyrotechnic products; 06 Sugar, sugar preparations and honey; 56 Fertilizers, manufactured; 52 Crude chemicals 
from coal, petroleum and gas. 


Ki 

iV’ii 

country 

IV’tl 

country 

IV’tl 

country 

\ipi\ 

country 



prod: 57 


prod:06 


prod:56 


prod:52 

1 

0.052 

USA 

0.216 

Mali 

0.332 

Brazil 

0.288 

Japan 

2 

0.044 

Tajikistan 

0.201 

Guinea 

0.304 

Bolivia 

0.279 

Rep. of Korea 

3 

0.042 

Kyrgyzstan 

0.059 

USA 

0.274 

Paraguay 

0.245 

China 

4 

0.022 

France 

0.023 

Germany 

0.031 

Argentina 

0.020 

Australia 

5 

0.021 

Mexico 

0.021 

Mexico 

0.017 

Uruguay 

0.013 

USA 

6 

0.018 

Italy 

0.021 

Canada 

0.009 

Chile 

0.012 

U Arab Em 

7 

0.018 

Canada 

0.018 

UK 

0.004 

Portugal 

0.010 

Canada 

8 

0.015 

Germany 

0.015 

Israel 

0.004 

Angola 

0.010 

Singapore 

9 

0.013 

U Arab Em 

0.015 

C d’Ivoire 

0.004 

Spain 

0.009 

Germany 

10 

0.012 

Qatar 

0.014 

Japan 

0.003 

France 

0.008 

New Zealand 


















